BIOMASS OF FIFTY CONIFER FORESTS AND NUTRIENT EXPORTS 
ASSOCIATED WITH THEIR HARVEST?! 


T. Weaver’ and F. Forcella’ 


ABSTRACT.— Biomasses of climax Rocky Mountain forests studied ranged from less than 50 to more than 300 
tons/ha. Total biomass was approximately 1.5 times the biomass of normally merchantable boles. When com- 
pared with conventional bole harvest, the nutrient exports associated with harvest of all aboveground parts in 
these stands would apparently be at least three times higher for nitrogen, six times higher for phosphorus, four 
times higher for potassium, and three times higher for calcium. 


The volume or biomass of a tree can be 
predicted from its diameter or from its di- 
ameter and height. Foresters routinely esti- 
mate the volume of logs using these dimen- 
sions and tables calculated on the 
assumption that tree trunks are truncated 
cones or parabolas (Forbes and Meyer 
1961). Weights can then be estimated by 
multiplying their volumes by the density of 
the wood. In most organisms the size of one 
organ system is directly related to the size 
of other organ systems (allometry, Sinnott 
1960), so one can predict the volumes or 
weights of whole organisms from simple 
measurements. This has been demonstrated 
for a wide variety of trees including those 
of the conifer genera Abies, Pinus, and 
Pseudotsuga (Ovington 1957, Ovington and 
Madgwick 1959, Kimura 1963, Baskerville 
1965, Whittaker and Woodwell 1968, Moir 
1972, Odegard 1974, Whittaker et al. 1974, 
and Long and Turner 1975). 

This paper presents a single equation 
which predicts total aboveground biomass 
of a tree, from its diameter, for any of the 
seven most common tree species of the 
northern Rocky Mountains. A comparison of 
this equation with others reported in the lit- 
erature emphasizes the similarity of diame- 
ter biomass relationships among conifers in 
general. Various shrubs may also have sim- 
ilar size-weight relationships (Weaver 1977). 
An equation predicting weight from mea- 


sures of both height and diameter is similar 
to published conifer D°H-weight relation- 
ships but less accurate than the diameter- 
weight relationship. Similar relationships are 
presented for weights of boles, branches, 
twigs, and leaves. 

The biomass-DBH relationship found was 
applied to the stand tables of Daubenmire 
and Daubenmire (1968) to estimate the 
aboveground arboreal standing crops of fifty 
forest stands representing the following cli- 
max associations: Pinus ponderosa-Festuca 
idahoensis, P. ponderosa-Symphoricarpos 
albus, P. ponderosa-Physocarpus malvaceus, 
Pseudotsuga menziesii-Symphoricarpos albus, 
P. menziesii-Physocarpus malvaceus, P. men- 
ziesii-Calaąamagrostis rubescens, and Abies 
lasiocarpa-Vaccinium scoparium. 

The nutrient content of the boles and of 
all aboveground parts was estimated so the 
nutrient export associated with conventional 
logging and complete tree utilization could 
be compared. 


METHODS 


Thirty-nine trees of six coniferous species 
and one broadleaf species were felled in 
southcentral Montana. The trees were sec- 
tioned into shoot material greater than 10 
cm diameter, shoots between one and 10 
cm basal diameter, shoots less than 1 cm 
basal diameter, and leaves. Each class was 
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weighed wet. Wet weights were converted 
to dry weights by multiplying by a con- 
version factor determined for each com- 
ponent of each tree. The conversion factor 
was computed as follows: dry weight (60 C 
to constant weight)/wet weight. The tree 
species and the numbers of trees per species 
sampled were Abies lasiocarpa (Hook.) 
Nutt., 7; Picea engelmannii (Parry) Engelm., 
3; Pinus contorta var latifolia Engelm., 4; 
Pinus ponderosa Dougl., 6; Pinus albicaulis 
Engelm., 7; Pseudotsuga menziesii (Mirb.) 
Franco., 5; and Populus tremuloides Michx., 
7. The trees sampled ranged from 4 to 54 
cm in diameter at breast height (D at 1.35 
m) and from 3 to 20 m in height (H). The 
trees came from a wide variety of environ- 
mental conditions ranging from dry to moist 
sites and from open to closed stands. A list 
of the trees and their habitat types, sizes, 
weights, and ages is available from the au- 
thors. 

After a preliminary graphical analysis, all 
data were pooled, and simple and multiple 
regression analyses were made to determine 
which models best predicted tree biomass, 
independent of species. These included lin- 
ear, curvilinear, and log linear regressions of 
12 dependent variables (including weight 
and logarithms of weights of total shoot and 
its components) against the independent 
variables D, D’, DH, and D°H. 

D°-biomass relationships were used to es- 
timate the aboveground arboreal biomasses 
of 50 forest stands for which Daubenmire 
and Daubenmire (1968) have published 
stand tables. Since their tables list tree num- 
bers by 10-cm-diameter size classes, tree 
weights were estimated by two different 
methods. 1) The weight of all trees in a 
plot were summed. These weights were esti- 
mated by assuming that the diameter of 
each tree was equal to that of the midpoint 
of the 10-cm-diameter class in which it was 
recorded. 2) The weights of all trees in a 
plot were summed. These weights were esti- 
mated by pooling all plot data for a given 
habitat type on the assumption that it rep- 
resented a single (abstract) forest; assigning 
trees in this sample to 2.5-cm-size classes 
from a balanced distribution of tree size vs 
tree number (Meyer 1953); and estimating 
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the weights of the trees as if their diameters 
were equal to the midpoint of the size 
classes. Regressions through logarithmically 
transformed data are said to underestimate 
biomass (Baskerville 1972); we used uncor- 
rected values because graphical analysis sug- 
gests that the uncorrected values represent 
our data better than corrected values do (cf 
Whittaker et al. 1974). 


RESULTS AND CONCLUSIONS 


Linear regressions of log,)D’ and 
log,,D°H against log weight generally pre- 
dict biomasses of Rocky Mountain trees 
well (Table 1 and Fig. la). (1) Both inde- 
pendent variables predict total biomass well 
(r? = 0.98 and 0.96, respectively). 
Log, )D°H may be the safer predictor if our 
regressions are used in habitat types where 
taller trees than ours may be found. (2) 
Log,)D°H vs log;y weight is the best pre- 
dictor for woody material with diameters 
greater than 1 cm or greater than 10 cm (r° 
= 0.98 and 0.95, respectively). This is 
likely so because the conic or parabolic 
models of stem volume are good, weights 
are strongly correlated with volumes, most 
of the weight is in the larger classes, and 
branch volume is both correlated with and 
overridden by stem volume. (3) Log, 9D’ vs 
log,) weight is the best predictor for leaves 
and branches with diameters less than ] cm 
(r? = 0.84 and 0.79, respectively). This may 
be so because the capacities of current 
transporting tissues is highly correlated with 
diameter (cf Shinozaki et al. 1964). A 
graphical analysis of the leaf weight-diame- 
ter relationship suggests that the relatively 
low r is partly due to between-species dif- 
ferences. (4) A regression of untransformed 
weights against D + D? was our best pre- 
dictor of branch weights (1 to 10 cm, r = 
0.97). This curve, not shown, is J-shaped be- 
cause in small trees the bole is in the 
branch-size class. (5) Regressions of untrans- 
formed variables (D, D’, DH, D°H) vs un- 
transformed biomass were not linear, were 
not used, and are not presented. 

Our regressions of log,;y)D’ or log,)D°H 
against log,, weight are remarkably similar 
to, but not identical with, those determined 
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in other regions (Fig. lb); predictions of 
weights from such relationships apparently 
become poorer as one considers species and 
environments increasingly different from 
those for which the relationships were de- 
veloped. The differences may be due to dif- 
ferences in form factor, wood density, 
branching patterns, and height-diameter 
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relationships. For example, the tendency of 
deciduous trees (Monk 1970, Whittaker and 
Woodwell 1968) to outweigh conifers of 
similar diameters could be due to the great- 
er wood densities and branchiness of the de- 
ciduous trees. 

The aboveground arboreal biomasses of 
50 stands are summarized in Fig. 2; the val- 
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Fig. 1. Biomass-dimension relationships for coniferous species. (A) Diameter-biomass relations for seven species 
of the northern Rocky Mountains: a = total biomass, b = stem biomass (diameters greater than 1 cm), c = bole 
biomass (diameters greater than 10 cm), d = twig biomass (diameters less than 1 cm), e = branch biomass (di- 
ameters 1-10 cm), and f = leaf biomass. (B) Published DBH-biomass and D*H-biomass relationships: a = Picea 
rubens (Whittaker et al. 1974), b = Pinus rigida (Whittaker and Woodwell 1968), c = seven species (Fig. la), d 
= Pinus contorta (Odegard 1974), e = Picea rubens (Whittaker et al. 1974), f = Pseudotsuga menziesii (Long 


and Turner 1975), g = Pinus contorta (Moir 1972), h 


Pinus sylvestris (Ovington and Madgwick 1959), i = 


Pinus rigida (Whittaker and Woodwell 1968), j = seven species (Fig. la), k = Abies balsamea (Baskerville 1965), 
= Abies veitchii (Kimura 1963), and m = Pinus sylvestris (Ovington 1957). 
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ues presented were calculated by assuming 
that tree diameters were equal to the mid- 
points of their 10-cm-size classes, i.e., by 
method 1. One may conclude (1) that total 
arboreal biomasses of climax conifer stands 
of the northern Rocky Mountains are gener- 
ally in the 50-300 tons/ha range and (2) 
that those of specific forest series or associ- 
ations generally have a narrower range: 
Daubenmire’s Pinus ponderosa-Festuca ida- 
hoensis stands had biomasses of 50-150 
tons/ha; his P. ponderosa-Physocarpus mal- 
vaceus, Pseudotsuga menziesit-Sym- 
phoricarpos albus, and Abies lasiocarpa-Vac- 
cinium scoparium stands had biomasses of 
100-250 tons/ha and his P. mensziesii-Phys- 
ocarpus malvaceus and P. menziesii-Calama- 
grostis rubenscens stands had biomasses of 
100-350 tons/ha. One metric ton/ha = 
1000 kg/ha. These estimates of standing 
crop do not include understory biomasses 
(perhaps 0.1 to 3 percent or possibly 7 per- 
cent of the arboreal biomasses: e.g., Oving- 
ton 1962, Whittaker et al. 1974, Tappeiner 
and John 1973) or belowground biomasses 
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(perhaps 20-30 percent or possibly 7-50 per- 
cent of the arboreal biomasses: e.g., Oving- 
ton 1962, Rodin and Bazilevich 1965, 
Whittaker and Woodwell 1968). 

Despite the exclusion of understory and 
underground biomasses, the standing crops 
reported are larger than most of those re- 
ported in the literature (30-70 tons/ha; Ro- 
din and Bazilevich 1965, Art and Marks 
1971, Whittaker et al. 1974). Reasons for 
the large standing crops observed are prob- 
ably two: Daubenmire selected stands much 
older than those commonly studied, and he 
selected the parts of those stands with the 
biggest trees. 

The accuracy of volume and biomass esti- 
mates can usually be increased slightly by 
recording narrower diameter classes (Meyer 
1953). To achieve this we pooled the data 
from all stands occupying a given habitat 
type on the assumption that it represented 
multiple stands of one (abstract) forest; de- 
termined balanced size distributions (Meyer 
1953), allocated trees to 2.5-cm-diameter 
classes, and estimated their weights by using 


Taste l. Regression constants and indices of error used in relating sizes and weights of important trees of the 


northern Rocky Mountains. 











Total Wood Wood Wood Wood 
Notes’ Tree >] cm >10 cm 1-10 cm <l] cm Leaves 
D? a 2 -1.05 -1.18 -1.40 0.60 -1.87 -1.91 
b 2 1.19 eet 1.25 0.75 1.16 1.13 
r? 5 0.98 0.97 0.89 0.81 0.79 0.84 
E 6 1.26 1.34 1.47 1.60 ES 1.91 
D-H a 3 Row —1.53 -1.84 20.73 -1.99 -2.08 
b 3 0.93 0.96 1.01 0.57 0.86 0.84 
t 5 0.96 0.98 0.95 0.76 0.69 ow 
E 6 1.36 Teor 1.30 les 2.64 2.19 
D+D a 4 =A +0.08 -29.36 + 29.37 Sale + 1.16 
b 4 oe -2.96 +0.71 -3.69 =0.57 + 0.39 
c 4 +0.85 +00.7] +0.49 +0.21 +0.15 +0.06 
i 5 0.96 0.96 0.94 0.97 0.63 0.56 
e 6 0.29 0.38 0.25 1.14 0.28 1.00 


Identification of symbols to the left. 


'LOg o WT = a + b logio D’ with weight in kg and diameter in cm. 





Logo WT = a + b log,, DH with weight in kg, diameter in cm, and height in m. 


WT = a+ bD + c D with weight in kg and diameter in cm. 


r is the correlation coefficient, r is the fraction of the variation explained by the regression. 


E is the antilog of the standard error of a logarithmic regression, expected errors lie between the predicted value x E and the predicted value/E. 
An estimate of relative error (e) is the standard error of a linear regression divided hy the average weight (y) used (Whittaker and Woodwell 1968). 
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the midpoints of those classes as diameters 
(method 2). Forest weights estimated by 
method 2 ranged from 80 to 98 percent of 
the mean of stand weights determined with 
10-cm-diameter classes (i.e., by method 1). 
They were: Pinus ponderosa-Festuca ida- 
hoensis 75 tons/ha, 80 percent; P. ponde- 
rosa-Symphoricarpos albus 155 tons/ha, 82 
percent; P. ponderosa-Physocarpos mal- 
vaceus 140 tons/ha, 98 percent; Pseudot- 
suga menziesii-Symphoricarpos albus 180 
tons/ha, 98 percent, P. menziesii-Calama- 
grostis rubescens 336 tons/ha, 91 percent; 
and Abies lasiocarpa-Vaccinium scoparium 
178 tons/ha, 93 percent. 
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Fig. 2. Aboveground arboreal standing crops of 50 
forest stands representing 7 habitat types: Pinus pon- 
derosa-Festuca idahoensis (A), P. ponderosa-Sym- 
phoricarpos alba (B), P. ponderosa-Physocarpus mal- 
vaceus (C), Pseudotsuga menziesii-Symphoricarpos alba 
(D), P. menziesii-Physoearpus malvaceus (E), P. men- 
ziesii-Calamatrostis rubescens (F), and Abies lasiocarpa- 
Vaecinium scoparium (G). Calculations include all 
trees with DBH’s greater than 5 cm; basal areas of 
stands 66, 101, 149, 173, 174, and 179 were recalcu- 
lated from stand table data. 
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Biomasses in the leaf, twig, branch, and 
bole compartments are very nearly (SE + 1 
percent) 9 percent, 13 percent, 13 percent, 
and 66 percent of the arboreal biomass in 
the mature forests studied. Total harvest in 
most Rocky Mountain forests should there- 
fore equal approximately 1.5 times the nor- 
mal (10 cm +) stem harvest (Fig. 2). 
Whittaker et al. (1974) give similar values. 
This is a consequence of both the allometric 
equations used and of the size class distribu- 
tions observed in the stands studied. 


Mineral nutrient export associated with 
total harvest is greater than 1.5 times that 
associated with normal bole harvest because 
needles and bark (which comprise a larger 
fraction of twigs and branches than of 
boles) are relatively rich in mineral ele- 
ments. Table 2 summarizes published values 
for nitrogen, phosphorus, potassium, and 
calcium contents for needles, twigs, 
branches, and boles. 


One can estimate the amount of a nutri- 
ent element exported from a harvested 
stand by adding up nutrient export associ- 
ated with each component exported (Table 
2). The nutrient export associated with har- 
vest of anv component equals the total ex- 
port x the proportion (percent) of the total 
due to the component considered (e.g., nee- 
dles or bole) x the nutrient content of that 
component. Since needles, twigs, branches, 
and boles comprise relatively constant pro- 
portions of total aboveground arboreal bio- 
mass (9, 13, 13, and 66 percent, respective- 
lv), we may conclude that nutrient export 
associated with total harvest—as opposed to 
bole harvest alone—is three to six times 
greater for nitrogen, six to nine times great- 
er for phosphorus, four times greater for 
potassium, and three times greater for cal- 
cium. Thus, total harvest of a Rocky Moun- 
tain conifer stand with a relatively low 
aboveground arboreal biomass of 100 
tons/ha might remove 160-460 kg/ha of ni- 
trogen, 20-40 kg/ha of phosphorus, 80-180 
kg/ha of potassium, and 120-500 kg/ha of 
calcium. Bole harvest alone would remove 
only 25-160 kg/ha of nitrogen, 2-6 kg/ha of 
phosphorus, 20-40 kg/ha of potassium, and 
40-200 kg/ha of calcium. Exports from for- 
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Tase 2. Typical nutrient contents of conifers.’ 
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Nitrogen Phosphorus Potassium Calcium 

kg/ kg/ kg/ kg/ 
% 100t % 100t % 100t % 100t 

Needles max 1.34 12) 0.22 20 0.78 70 1.05 95 
min 0.84 76 0.08 1 0.28 25 0.15 14 

Twigs max 0.60 78 0.07 9 0.30 39 1.10 143 
min 0.31 40 0.05 a 0.14 18 0.27 35 

Branches max 0.77 100 0.06 8 0.25 32 0.54 70 
min 0.15 20 0.02 3 0.12 16 0.28 36 

Bole max 0.24 158 0.01 7 0.06 40 0.30 198 
min 0.04 26 0.00 2 0.03 20 0.06 40 

Total max 457 43 181 506 
min 162 18 79 125 





'Percentage composilions reported are the highest and lowest values reported for a variety of conifers (Larix, Picea, Pinus, and Pseudotsuga) in papers 
by Odegard (1974), Overton et al. (1974), Stark (1973), and Rodin and Bazilevich (1965). 


*Kg/100t contents refer to 100 tons of aboveground slanding crop. They were calculated by multiplying 100 tons of standing crop x the approximate 
percentage of the standing crop contributed by the organ (9 percent needles, 13 percent twigs, 13 percent branches, and 66 percent boles) x the per- 


centage of the organ contributed by the nutrient considered (N, P, K, or Ca). 


est twice to three times as productive 
would be twice to three times as great. 

One wonders whether the forest product 
value of ‘slash’ will exceed its fertilizer val- 
ue. If inputs from nitrogen fixation process- 
es or from decomposition of rocky sub- 
strates are small, the fertilizer value of 
branch material with diameters less than 10 
cm might exceed its forest product value. 
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